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g . Abstract 

FT' 

Magic-angle spinning (MAS) solid state nuclear magnetic resonance (NMR) 

"X' 

' spectroscopy is shown to be a promising technique for implementing quantum 

computing. The theory underlying the principles of quantum computing with 
nuclear spin systems undergoing MAS is formulated in the framework of for- 
malized quantum Floquet theory. The procedures for realizing state labeling, 
state transformation and coherence selection in Floquet space are given. It 
suggests that by this method, the largest number of qubits can easily surpass 
that achievable with other techniques. Unlike other modalities proposed for 
quantum computing, this method enables one to adjust the dimension of the 
working state space, meaning the number of qubits can be readily varied. The 
universality of quantum computing in Floquet space with solid state NMR 



is discussed and a demonstrative experimental implementation of Grover's 
search is given. 
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I. INTRODUCTION 



As early as the late 1950's, Landauer and Bennet et a/.^'^'^'^'^'^, investigated the ef- 
fects of physical laws on computing, such as the reversibility of a computing operation and 
the minimal energy required to transmit a bit of information. Feynman^, on the other 
hand, was studying the fundamental limitations of quantum mechanics on the capacity of 
(classical) computers. The most important question in these works was what would it hap- 
pen if computing logic is not presumably given but rather determined by physical laws, 
particularly, quantum mechanical laws? With the rapid development of very large scale 
integrated circuitry technology, above question seemed to become important in the early 
1980's; that can be rephrased as, what would it happen if the chip size were made so small 
that one chip contains very few, even just one impurity electron. That background of sci- 
entific development initiated quantum computing research. However, quantum computing 
was basically dormant in the decade of the 1980's. It has since gained increasing attention 
once the power of a hypothetical quantum computer was revealed, particularly, through the 
works of Deutsch^'^'^°'^^'^^ Shor^^'-*^^, Lloyd^^ etc. Deutsch^'^'-^^ showed that genuine and mas- 
sive parallelism can be achieved. Lloyd^^ proposed a quantum computing prototype that 
has subsequently been followed. Shor^^ demonstrated the power of quantum computer in 
solving the famous and all-important problem in number theory and public key cryptog- 
raphy system, i.e. the prime factoring of large integers. Shor et a/.^'''^^'^^, Gottesman^®'^^, 
Steane^'^'^^'^^, Schumacher^^and PreskilP^ and others invented a variety of quantum error 
correction schemes that are crucial to the realization of long-time quantum computing. 
Since then, theoretical publications have appeared with increasing frequency, encompassing 
almost every aspect of computing theory(for review, see, e.g.,^"^'^^). Remarkable progress in 
experimental implementation and model proposals also has been made in utilizing an exten- 
sive repertoire of sophisticated experimental techniques including atomic interferometry^^, 
quantum electrodynamic cavity^^'^^'^^'^^'^^, ion trap^^'^^, polarized photons^^, nuclear spins 
embedded in an electron system in the quantum Hall regime^®, quantum dots^^, Joseph- 
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son junction electrons in liquid helium , nuclear spins in doped silicon devices , single 
Cooper pair^^ Rydberg atom^^ and liquid NMR^^'46,47,48,49,50,5i,52,53,54,55,56,57,58,59,60,6i j^jf. 

fering from other techniques, the NMR prototype uses bulk samples hence an ensemble of 
nuclear spins rather than pure quantum mechanical systems. Among the above experimen- 
tal prototypes, NMR is certainly the most promising, to date: all above methods except 
NMR can only simulate a single quantum logic gate such as controlled NOT gate, but NMR 
can do much more than that, e.g., it can simulate a quantum network such as the per- 
forming of simple arithmetic operations, and a quantum computer that can execute simple 
quantum algorithms^®'^^'^^'^°'^^'^^'^^. The NMR method offers the first realizable quantum 
computer operating with more than two qubits, thus providing for the first time a quantum 
computer with error correction capacity^^. All these demonstrations used liquid state NMR 
spectroscopy because of its natural high resolution. While the progress has been remarkable, 
one severe difficulty with the NMR quantum computer is the exponential loss of the signal 
sensitivity with the increase of spin numbers (hence usable qubits) in the working molecule. 
It is clear that establishing a NMR quantum computer with a capacity of over ten qubits is 
rather challenging^^, if not impossible even though a host of sensitivity-enhancing techniques 
are available^^. 

In this paper, we present an alternative, probably more advantageous, method for performing 
NMR quantum computing, that is, quantum computing based on solid state NMR involv- 
ing rotating samples at an angle of 54.74°, the magic-angle to the apphed magnetic field. 
This so-called magic-angle spinning (MAS)^^'^"^ NMR can be well formulated using Floquet 
theory^^'^^'^^'^^'^^'''^'''^'''^'''^'^^. From the point of view of quantum computing, the Floquet 
description offers a method to augment the state space, almost infinitely. In practice, nev- 
ertheless, the size of the space is restricted by the signal sensitivity. However, as shown 
in our theoretical analysis elaborated below, this size can be easily made much larger than 
that realizable in liquid NMR studies. For quadrupolar nuclei, the sidebands produced by 
the rotating polycrystalline samples can be as many as thousands or even more'^^''^^, mean- 
ing usable qubits can be easily achieved, even in excess of 10 merely by using conventional 
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NMR techniques, although in this case, the manipulation of quantum states and coherences 
is more complicated than for spin-1/2 systems. The other obvious advantage of solid state 
NMR is that the number of spins in a sample is usually much larger than in a liquid sample 
of the same size, meaning a significant sensitivity gain. Our paper is arranged as follows. 
In the Section II, the theory underlying quantum computing with sohd state MAS NMR 
is described; this is the foundation of the subsequent sections of paper and future work. 
Particularly important are the definition of the pseudo pure state in Floquet space and its 
connection with quantum computing. While the theoretical framework applies to nuclear 
spin-1/2 systems as well as to quadrupolar nuclear spins, the remainder of the paper will 
focus on spin-1/2 systems with chemical shift interactions. Section III estabhshes the corre- 
spondence between a pseudo pure state and its spectral representation. This is essential to 
the read-out function in NMR quantum computing because the directly detectable signals in 
NMR arise from the single quantum coherences. The theoretical derivation of the spectral 
signal is demonstrated in Appendix A. The preparation of pseudo pure state is crucial to 
quantum computing and this is discussed in Section IV. Three different methods are con- 
sidered. Section IV analyzes the universahty of MAS solid state NMR quantum computing. 
In Section VI is demonstrated the implementation of an important quantum computing 
algorithm, namely, Grover's search, on a solid state NMR quantum computer. The major 
points of this paper are summarized in the final Section. 



II. QUANTUM COMPUTING IN FLOQUET SPACE 

A. Quantum Floquet theory of solid state NMR 

A periodic time-dependent Hamiltonian such as that for a nuclear spin system in a 
polycrystalline sample undergoing rotation at the magic-angle in a static magnetic field is 
best described employing Floquet theory. Here we summarize the well-developed theory from 
the perspective of MAS NMR and its significance to quantum computing. Most generally. 
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the evolution of the density matrix, p{t), of a spin system can be written as 

p{t) = fe-^/o <^*'^(*')p(0)e^/o '^*'^(*'). (2.1) 

It follows, therefore, that evaluating the evolution operator U{t) — e~^->o^^^^^^ is a central 
part of spectral lineshape calculations. The straight-forward procedure is to use the multi- 
step method which divides the time interval, (0, tc), where tc is the period of the Hamiltonian, 
into N equal steps and then one calculates each step by approximating its Hamiltonian as 
being time-independent 

where t — ntc/N . This usually involves the diagonalization of each instantaneous Hamilto- 
nian H{ti). Floquet formalism^^'^^'^^'™'^-^'^^'^^'^^, on the other hand, focuses the calculation 
of the evolution operator on computing the Floquet Hamiltonian Hp by introducing Floquet 
states \rn > where r is the state index of H and n is the mode index^''''^^ 

< pm\HF\qn >= + noUcSpqSmn (2.3) 
where hp^ are the Fourier components of the time-dependent Hamiltonian 

^^P. = E^V^* (2-4) 

k 

The evolution operator then can be calculated from the following expression 



n=—oo 

oo 

= E E < P^l^r >< > e-'^'^^-'^'^^K (2.5) 

f n,fc=— oo 

where the index r runs over the Hilbert space defined by H . For the sake of generality, the 
Hamiltonian is assumed to be anisotropic, i.e. H = H{a, (3,j,t) where a, (3 and 7 are the 
Euler angles describing the interaction tensors relative to the laboratory frame (they can be 
set to zero for solution NMR cases). |A" > and A" are the eigenstate and eigenvalues of the 
Floquet Hamiltonian, respectively 

6 



Hf\K > = K\K > (2-6) 

= A" - nujc (2.7) 

Eq.(5) can be rewritten as^^'^^'^^ 

U{t) = P(t)e-*^*P(0)-^ (2.8) 

where Q is traceless and diagonal with diagonal elements qr, and P{t) the Floquet amplitude 
is defined as 

Pp,=<Pn\X'^> (2.9) 

At first sight, the estimated magnitude of U{t) can be made almost exact because the values 
of t can be chosen in arbitrarily small increments. However, in most realistic cases, the 
Floquet Hamiltonian Hp cannot be solved exactly: it requires the use of a perturbation 
expansion, or equivalently, matrix diagonalization. As has been shown, the order of the 
expansion series, or the chosen dimension of Hp presents a bound to the accuracy of U{t). 

If we further define density matrices, observable operators and evolution operators in the 
Floquet basis, a formalized Floquet theory^^ can be formulated. Specifically, we define 

U''{t) = ^n,mUn-m{t)\n >< m\e-'^^' (2.10) 

where Un{t) are given by 

C/(t,to) = S„[/„(t)e-^"'^'-*° (2.11) 
Then the density matrix can be found as 

a{t) = E„,^ < n\a''{t)\m > e-("— (2.12) 

where satisfies the Liouville equation 

^-t[H',a'{t)] (2.13) 



dt 



and the initial density matrix at time i = is given by 

< n\a^\m >= 5„,^(7(0) (2.14) 

The observable operator is defined as 

= En^^An_m\n >< m\ = E„,„A„|n >< n - m\ (2.15) 

where are the Fourier components of A[t): 

A{t) = E„^„e-*"'^'-* (2.16) 

Note that definition of the Floquet Hamiltonian in the original work^^ is different from 
Eq.(3). The formalized form of the detection observable D then is easily found from the 
requirement 

S{t) = Tr[Da{t)] = Tr[D^a^{t)] (2.17) 

to be 

= j:n,mD\n >< m\e-'^">'^^ = S„,^i^|n >< n + m\e-''^'^^ (2.18) 

Because of its unified and has a concise form, the formalized Floquet theory will be used 
throughout the work. 

B. Floquet pseudo pure state 

With the above formahzed Floquet theory, the pseudo pure state or effective pseudo pure 
state introduced by Cory el al!^^ can be extended to Floquet space as follows: 

where n is the number of spin-1/2 nuclei, k is the (effective) dimension of the "mode" space 
and I is the identity spinor whose matrix form is an nk x nk identity square matrix. It is 
easy to verify that the above definition of the pseudo pure state satisfies the three criteria 

8 



given by Cory et al.^: i.e., "if^ defines a pure density matrix and vice versa; evolves 
according to the same unitary transform governing the evolution of a pure density matrix and 
the measurement value of an observable operator over and that of the same observable 
operator over a pure density matrix differ in a trivial constant only. Therefore, the pseudo 
pure state in Floquet space can be used to "emulate" quantum computing. 

To specify the Floquet space, in the following sections we will focus on the chemical shift 
interaction. The Euler angles system is defined as follows: the principal axis direction of the 
chemical shift tensor in the rotor systems is determined with (a, /?, 7) while the rotor system 
is specified by {cUrt, d-i 0) relative to the laboratory frame, where ujr is the sample spinning 
speed. The most interesting case in solid state NMR as well in this work is when the sample 
spins at the "magic angle", i.e. the spinning axis is tilted P — Prn — 54.74° with respect 
to the static magnetic field, which is called magic-angle-spinning(MAS). The chemical shift 
interaction Hamiltonian then can be written as^^'^^ 

Hcs = -Wo + 5P2(cos^)[P2(cos/3) - |sin2/3cos27] + ^^^(t)} (2.20) 

where 5o is the isotropic chemical shift plus the RF offset, 5, rj are the anisotropy and 
asymmetry parameters of the chemical shift tensor, respectively. Denote the three principal 
values of the chemical shift tensor as au, C22, o"33 and assume the convention an > (J22 > (733, 
then there are relations Sq — |((7ii+(T22+(J33), S — ctq— (T33, 77 = (cii— C22)/(^- It is noteworthy 
that at the magic angle, all the anisotropic terms in Eq.(20) disappear, meaning registered 
spectra are free from line broadening caused by chemical shift interactions. This is the most 
important principle in high resolution solid state NMR. The time-dependent term ^(t) in 
Eq.(20) is given by 

^(t) = CiCOs(u;rt) + Sisin(ujrt) + C2COs{2ujrt) + S2S\n{2ujrt) (2.21) 

where 

Ci = ■^sin2^sin/3[— cos/3(77COs27 + 3)cosq; + 77sin27sinQ;] (2.22) 
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Si — -sin2^sin/3[cos/3(r7COs27 + 3)sinQ; + ?7sin27cosQ;] 

13 77 

C2 = -sm^9{[-sm^(3 — -cos27(l + cos^/9)]cos2q; + 77COS/9sm27sm2Q;} 

1 3 7? 

'S'2 = -sin^^{[— -sin^/? — -cos27(l + cos^/?)]sin2Q; + r7cos/?sin27cos2Q;} 

Zi Zi Zi 

Prom Eqs. (20,21), the Floquet Hamiltonian Rq^ can readily be found from Eqs.(3,4) 
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<^n\HQg\^n> — nujr 5uj < ln\Hl;g\ln >— nujr -\ — 5a; (2.23) 

/3 /3 
< On|//^s|On± 1 > = -^(^csCi < ln|//^5|ln± 1 >= ^fcsCi 



< Qn\H^s\Qn ± 2 > = -^^csCs < ln|i/^5|ln ± 2 >= ^(5c5C2 



All the other elements are zero. The chemical shift interaction is a typical "inhomogeneous" 
interaction, i.e., it's Hamiltonian at different times is always commutable rendering it un- 
necessary to perform the time ordering operation in the calculation of the unitary evohition 
operator. This is a very important property that helps one analyse the evolution of (pseudo 
pure) quantum states and simplify the design of quantum computing gates. Other important 
inhomogeneous interactions include the electric quadrupolar interaction and heteronuclear 
dipolar coupling which will be discussed in future work. 

Prom above paragraphs, some important implications of the applications of Ploquet 
space and solid state NMR to quantum computing in Floquet formalism are summerized as 
follows: Pirst, Ploquet space is dimensionally adjustable, i.e., changing the sample spinning 
speed ujr can augment or reduce the effective dimension of the space meaningful for quantum 
computing. Moreover, contrary to a usual quantum mechanical system, its dimension is not 
necessarily a power of 2. Second, in solid state NMR quantum computing, the Hamiltonian 
hence the "operation" can be controlled both with the RP field and sample spinning speed, 
which provides more flexibility than solution NMR quantum computing. Third, opposite to 
solution NMR quantum computing which is not satisfactory at low temperatures (below the 
melting point of the sample used) , sohd state NMR is usually more sensitive at temperatures 
as low as possible. 



10 



III. SPECTRAL REPRESENTATION OF FLOQUET STATE 

A. Signal readout 

One of the most essential functions of computing is that the output can be read out. 
In this section, we give an operational procedure on how to "read out" a pseudo pure 
state of a solid state NMR quantum computer. In hne with standard NMR spectroscopy 
which measures the induction voltage caused by the transverse magnetization vector of the 
spin ensemble, this is done by observing the spin ensemble (which is in a pseudo pure state). 
There are numerous " read" pulses available, but for simplicity and without losing generality, 
we use single 90° pulse in this work (This is sufficient for chemical shift interaction but 
for quadrupolar interaction an effective observation may demand more complicated pulses, 
which will be discussed in the future). 

The energy levels of a spin- 1/2 system in spinning sohd NMR are labeled as shown in 
Fig. 1(a). The first index is spin angular momentum quantum number and the second one 
the mode. The readout function of an output state is given in Fig. 1(b). Therefore, the FID 
(free induction decay) signal of a pseudo pure state \pm > can be given as 

Snt) = Tr[I^U^sU[o'^'^'-Ui,-'u^s~'] (3-1) 

where U^g and UgQ are the Floquet evolution operators corresponding to chemical shift 

interaction and the 90° pulse, respectively, a^^"^ = \pm >< pm\ is a pure state in Floquet 
space. is the observable operator defined by 

la = ^m,n\n > la < n + m\e''^^^ a^x,y, ±, etc. (3.2) 



B. The expressions for Ug^ and U^g 

We assume the 90° pulse is along the -x direction in the rotating frame, the RF Hamil- 
tonian is written as 
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Hrf — LOiIx 



(3.3) 



where Ui is the RF field strength and is the Pauh matrix. Above equation leads to the 
following Floquet Hamiltonian 



H^j: = ct;iS„|n > Ix < n\ + nujrl\n >< n 
^2^^ ® ® ■■■'2 



(3.4) 



n\ 



The general expression for the evolution operator of the RF interaction is then given by 



U^j^(t ) = Q~^^rf^P — Q-i[i^i'^n\n>Ix<n\+nuirl\n><n\]tp 

where tp is the pulse width. In explicit matrix form, Eq.(28) is 



(3.5) 



Ix 



(3.6) 



V 



If the effective dimension of the mode space is K and the condition Kuj^tp — > is satisfied, 
above equation is reduced to 
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u, 



Ix 



(3.7) 



V ■/ 

which is a useful simphfied expression. The exphcit expression of the evolution operator 
of chemical shift interaction Hamiltonian can be found from Eqs. (20-22). Specifically, for a 
spin-1/2 system, from the Hamiltonian Eq.(20), we have 



TT (+ 4- \ -iJ*dt{u>cs+So)Iz 



(3.8) 



— ^n^n 







i[^(t-3to)+%^] 



e-L 2 



I 



Where the expansion coefficients A„ can be found to be = |F„|^ with^^ 



ZTT JO 

Comparing Eq.(31) with Eq.(ll), we have 
which, when = is chosen, is reduced to 



■^sinoi>— -^008(21)+ sin2{T^— cos2(il 



(3.9) 







gi[^(t-3to)+^] 



(3.10) 



6 2 







.So+nujr. 



(3.11) 



which can, in terms of its matrix elements, be denoted in a more concise form as follows 
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U'cUt) - Ar.e-'''^^^ p = 0,l (3.12) 

where p — 0(1) corresponds to spin state |0 >= | + | > (|1 >= | — | >) and €p — p — Sp^ 
with 6p^o the Kronecker function. 

C. The effects of Uqq and U^g on pseudo pure state 

Suppose we have prepared a pseudo pure state \pm >by state labehng techniques'^ '^^'^^ 
(detailed in Section IV for sohd state NMR quantum computing). The effect of the 90° RF 
pulse on the state is 

\pm > E„,,e-^"'^'-*^|n > [7^°^ < l\pm > (3.13) 

90° 

Because the only non-zero component of is Uq^^ and notice that |0 > — ^ (|0 > > 
)/\/2, |1 >-^ (|1 > -|-i|0 >)/\/2, the above equation is simply 

90° 1 

\pm > — ^ — > (|p > +i\p — >) (3-14) 
v2 

The effect of the chemical shift interaction on a pseudo pure state can be found as follows 

\pn >^ U^slpri > = ^k,ie-"''''-'Ucsk-i\k >< l\pn > (3.15) 
= i:ke-"''^^'Ucsk-n\pk > 

The combined effects of the RF pulse and the chemical shift interaction are therefore given 

as 

UEsUL\wi >= -^EkC-^'^'^'-'Ucsk-milpk > +i\p - epk >) (3.16) 
Eqs. (37-39) are the basic equations important to the calculations in the following sections. 
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D. The spectral representation of the readout signal 



The signal Eq.(24) can be decomposed into two terms S^{t) = S^{t)+iS^{t) where 
Sl"^{t),S'!y^{t) are obtained from Eq.(24) by replacing /+ with and ly, respectively. We 
will give the results here (the derivation of S^"^{t) is shown in Appendix A and it equally 
applies to the calculation of S^{t)). 

S^{t) = i-EfcAfc_^ep[e-^'^(*^-™)'"'-* - e^^^^^"'")^'^*] (3.17) 

whose spectrum is obtained by Fourier transformation 

/^"*(a;) = i^EkAk-m^p[S{uj - ep{k - m)ur) - 8{u + ep{k - m)ur) (3.18) 

For a given system, the largest mode number is fixed, say, which is restricted by the 
sensitivity limit. Then the above equation means that the sideband manifold consists of the 
following bands: [—K — m, —K — m+1, X — m]ep and [—K + m, —K+m + 1, K + m]ep. 
The amplitude of each band depends on the value of p: for different p, there is a 180° phase 
factor difference. Therefore, a unique one-to-one correspondence is established between a 
pseudo pure state and a spectral representation. The signal <S'^™(t) (or /^™(a;)) contains two 
groups of bands with opposite signs in intensity. If quadrature detection is used, the signal 
is found to be 

^r(^) = ^kAk-mepe-''"^'-'^^'''^' (3.19) 

and its spectrum is readily found to be 

/^'"(^) = EkAk-mepS{uj - ep{k - m)w,) (3.20) 

which contains only one group of sidebands given by indices [—K + m, —K + m + 1, T^' + 
m]ep. As a demonstration, we choose the six-level system shown in Fig. 2 in which \K\ — 1. 
Then the spectral representations of all the pseudo pure states are given in Fig.3 and Fig.4 
for single crystal and powder samples, respectively. The one-to-one correspondence between 
state and spectrum is obvious. The pseudo states are differentiated by amplitude distribution 
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or by a phase factor or both. The spectra are absorptive for powder sample but for single 
crystal, they may be dispersive but in the later case, the relative phase differences can still 
be used to unambiguously identify the particular pseudo pure state of the solid state NMR 
quantum computer before the readout RF pulse was applied. 

IV. STATE LABELING 
A. General 

State labeling is the unique feature of ensemble quantum computing because a pure state 
is not naturally available in an ensemble. To carry out quantum computation, one must first 
"purify" the ensemble so that it can be regarded as being in a pure state. The problem of 
state labeling in solution NMR has been well addressed. There are two types of labeling: 
one is temporal and another spatial. The former uses proper time averaging and the latter 
uses spatial averaging, of the density matrices, to construct a pseudo pure state. Here we 
show how these methods can be extended to Floquet space quantum computing. 

Given an initial density matrix, p^{to), state labeling renders it a pseudo pure state and 
a quantum computation task can be undertaken starting from it as follows 

P^C/V(to)C/^"' = C\<P^ X 00^1 = (4.1) 

where C is the computation operator and c is a constant coefficient and a certain 
preparation (not unitary in most cases). 

State labeling from the initial (thermal state) density matrix given by < n\a^ {0)1171 >— 
c(0)5o,r„5o,n is a trivial task because there are only two non-zero terms for spin-1/2 systems. 
Either with phase cycling (two experiments) or by applying a gradient field, a specified 
state can be chosen using the routine techniques. However, this is not the general case 
because, first, the initial density matrix in formalized Floquet theory can be chosen as a 
different form < n\a^{0)\m >= (j{0)6m,n and secondly, state labeling may need to start with 
a density matrix other than the thermal equilibrium form. In addition, certain computations 
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may require a pure state with the mode indices not equal to zero. Therefore, in the following 
two general techniques will be given which can be employed for any type of mixed state. 

B. Multi-pulse techniques 

Prom Section III, a pure state corresponds to certain peak profile in an MAS spectrum. 
The preparation of a pseudo pure state, therefore, is equivalent to constructing a subspcc- 
trum with specific peak profile. Over the past decades, there have been developed an array 
of techniques in sohd state NMR to manipulate the MAS spectral peaks, such as TOSS 
(total suppression of sidebands), PASS (phase adjusted spinning sidebands), IRS (isotropic 
rotation sequence) and their combinations and two-dimensional extensions etc^®. As an 
example, we demonstrate that the 2D-PASS sequence is a satisfactory technique for state 
labehng. The sequence consists of five tt pulses and is shown in Pig. 5. The parameter 
(pitch) is used to characterize the separations between each adjacent pair of tt pulses. In 2D- 
PASS experiments^, © also represents the first dimension and it is changed systematically 
according to prescribed separations between pulse pairs such as 

2S^=i(-l)''e^"*^« + 1 - (-i)^e^"^(®+^^) = (4.2) 
-2E^=i(-l)"+<^^, + ^T-0 

where n is the total number of the sidebands in the MAS spectrum. The meaning of other 
parameters is shown in Pig.5. Each sideband (or central band) can be extracted from the 
projection along the first dimension. When necessary, the intensity of each band can be 
modified by adding a pair of pulse %Q°_^dx at the end of the PASS sequence, where the value 
of is determined by the desired intensity of the band. The peak profile of a pseudo pure 
state can then be constructed by adjusting © and 9x systematically. Por example, given a 
peak profile Ak-,K — 1,2, ...,n where Ak is the intensity of K-th sidebands. The second 
dimension signal (PID) in a 2D-PASS experiment is given by 
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S(t2) = E^_^a(*^)e-^'=®e^(''°+'='^'-)*2 (4.3) 

where the exphcit expression of a^*^^ can be found in Ref.[77]. The simplest peak "basis" can 
be chosen as 

Se-5®(t2)x| = AKe'^^°+^''^^'' (4.4) 

where xk {— sinOx) are to be determined. By summing all possible values of © the above 
equation can be simplified as 

Eea^^^e^'^Qxl = S^xAk (4.5) 

Therefore, the values of xk can be found out by solving the above linear equations 

X = A-^A (4.6) 

where the matrices A = a('=)e**^® and A = S^kAk- It can be easily verified that A~^ exists. 

Therefore, any subspectral profile can be constructed with this sequence: 90° — [ASL] — 
90° — 9_x- The number of the values of © is determined by the number of the sidebands in 
the MAS spectrum and in principle there are no restrictions on the values of 0. Therefore, 
this technique can readily deal with a system with as many as hundreds or even thousands 
of sidebands. 

To reduce the number of steps to a practically acceptable level, the density matrix may 
be simplified before the standard state labeling techniques is invoked. We believe that the 
selective excitation techniques such as Dante pulse sequence^*^, can be incorporated in the 
preparation of the density matrix for state labeling. 

C. Gradient field selection 

Temporal labeling methods as discussed above are easy to implement and retain the 
sensitivity of the whole sample, but usually require a number of experiments. This re- 
duces computing efficiency and in the least-promising cases might conceal the advantage of 
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quantum parallelism^^'^^. Therefore, temporal labeling is best suitable for low-qubit imple- 
mentations. Another technique, i.e., spatial labeling was proposed by Cory et a/.^^'^^ which 
may or may not use phase cycling. Here this method is extended to the solid state MAS 
NMR case. 

In a magnetic gradient field, the Floquet energy levels are dispersed along the direction 
of gradient. The effect of a gradient field on the energy levels is the same as in normal 
Hilbert space. Without losing generahty, we consider the case where the gradient field is 
along the z-direction (parallel with the static magnetic field) . We use the " sandwich" pulse 
sequence ^^1(^01) — 90^ — ^22(^02) where tci,tG2 are the gradient pulse lengths and 90° is 
the RF pulse, as shown in Fig. 6. The evolution operator during the gradient pulse can be 
from Eq.(ll) 

UG{t,to)^e ■'*o (4.7) 



^izG4^it-3to)+"-f] 

S„C/Gn(i)e'""'^*° 



V 



from which the evolution operator in Floquet space can be easily found using Eq.(lO). 

With any given initial state p^{to), the final density matrix after the sandwich pulse 
sequence is given by 

p^{t) = U^,{h,tM,US,{tGuto)p^{^^^^ (4.8) 

where t2 — ^02 + ^o, = + tai +tp. It is easy to show that above equation only provides a 
constraint on the the Floquet space if the gradient field is time independent. Using Eq.(lO), 
we can find the condition for state labeling from Eq.(51) as 

epku}rGz2tG2 + SqlUrGzltGl = (4.9) 

where k, I are integers. The detailed derivation of above equations is given in Appendix B. 
By setting Gzi,Gz2,tGi and tG2 properly, the desired situation that only one term survives 
the pulse sequence can be realized, thus a pseudo pure state is prepared. 
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Obviously, the gradient field method can also be incorporated with selective excitation 
but its merit is unknown. In the following, we will focus on the techniques that do not need 
a gradient. The thorough investigation of the application of a gradient field in solid state 
NMR quantum computing will be presented elsewhere. 

V. UNIVERSALITY AND GATE DESIGN 

The above procedures for implementing state labeling can be extended to construct basic 
operating matrices, i.e., elementary quantum logic gates. The universality can be ensured if 
all possible operations can be realized with a set of elementary gates. The central problem is 
how to realize the operations that transform any state into any other state. As an example, 
we show that the 2D-PASS sequence can be used as a basic pulse block in constructing 
universal logic gates. Prom "peak manipulation" point of view, if all peak profiles can be 
realized starting from any peak profile, the gates thus constructed are universal. 

Here we give a general discussion on the relationship between peak manipulation and 
the construction of complete unitary transform (universal gates). Starting from the initial 
density matrix p^(0), two pseudo pure states can be prepared as follows 

pf = P,V(0) (5.1) 

= p,V(o) 

Thus the unitary transform that maps state p{ to state pf can be realized with 

U[, = P[Pr' (5.2) 

Suppose there are L Floquet states that are usable for quantum computing where L — 
Ld + Lu where Ld, I/„ are the number of states corresponding to spin quantum number 
|, — |, respectively. Its is easy to find from the energy levels that the total number of peaks 
from this energy level manifold is L — 1. Therefore, for a unitary matrix of M x M dimension. 
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M + 1 Floquet states are needed. The unitary matrices can then be determined from the 
peak manipulation matrices which constitute a complete set: 

UL = S,,a-"P;7" (5.3) 

The unitary matrices thus constructed necessarily form a complete set because peak manip- 
ulation matrices are complete. 

The experimental realization of the universal operation can be implemented with, e.g., 
2D-PASS sequence [ASL] shown in Fig. 5. For example, pulse sequence 

90^^ - 90° - [ASL] - 90° - e^Z\ where [ASL]-^ is the time-inverted version of the ASL 
sequence shown in Fig. 5, gives the pulse sequence that transforms a pseudo pure state 
\pm >< pm\ to another |gn >< qn\. With these transform matrices, therefore, any possible 
operations compatible with the system are realizable. 

This method, perhaps rather forceful, is workable in practice, at least for low qubit cases. 
However, we point out that there may exist more efficient pulse sequences that can realize 
above transformations and that is to be the major goal of our subsequent efforts. 

VI. EXAMPLE: GROVER'S SEARCH 

In this section, we demonstrate the use of the theoretical formalism developed in the 
preceding sections to implement experimentally a quantum algorithm. We will use Grover's 
search^^ as an example. 

Grover's search consists of four steps^^: (l)the preparation of pseudo pure state; (2) HW 
transform; (3) conditional flipping; (4) average about the mean; where step (3) and (4) are 
repeated 2y/N/7r times for an A^— item search. When using the 2D- PASS sequence, steps 
(1-3) can be incorporated in one experimental step with the initial density matrix as the 
input and an equi- amplitude superposed state as the output. In the output state, the phase 
difference between the flipped bit and the rest is 180°. The step (4) corresponds to the 
transform matrix with elements Uij = | — 5ij, each of which can be implemented according 
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to Eq.(55). Generally, the corresponding pulse sequence can then be designed as follows: 
90° - [ASL] - 90° - e^ll - ef^ - 90^^ - [ASL]-^ - 90° - [ASL] - 90° - - Acquisition with 
l,m,n = 1, 2, 3, 4... This means the number of experiments is / x m x n which would soon 
become impractically large. However, above construction of the pulse sequence may reduce 
the number of experiments substantially. For example, for the 2-qubit search experiment 
demonstrated in this work, however, only four experiments are required for each search. The 
Floquet energy level diagram and four "working" states are illustrated in Fig. 7 (A) and 
(B), respectively. 

The experiments were performed on a Bruker MSL-200 NMR spectrometer. The ^^C 
resonance frequency was 50.3 MHz. The sample was hexamethylbenzene (HMB) which was 
spun at the magic-angle with a spinning speed of 5 kHz. The aromatic carbon atom which 
has a chemical shift anisotropy of about 100 ppm is the spin we used for quantum computing. 
The methyl carbon has a small chemical shift anisotropy and its peak is used as the phase 
standard of the spectra. With the spinning speed used, there are two sidebands that are 
clearly observable for the aromatic carbon. Based on pulse sequences and the procedures 
discussed in the previous section, four search results were obtained and are shown in Fig. 8. 
As seen from Fig. 7(B) and Fig. 8, the theoretical prediction is in good agreement with the 
experimental result. 

VII. CONCLUSIONS 

In this work, magic-angle spinning sohd state NMR QC theory is developed based on 
formalized Floquet theory. Because the mode space is controlled by the sample spinning 
speed, the realizable number of qubits is changeable and can be made as large as the task 
demands (the ultimate limitation comes from signal to noise ratio) . The techniques required 
for state manipulations (labeling, coherence selection etc) are analyzed. The spectral repre- 
sentations of state is given that is crucial to readout of QC registers. The basic QC gates 
are demonstrated and an important QC algorithm is shown to be realizable in solid state 
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quantum computer. Based on these theoretical analyses, it has become clear that solid state 
NMR may become a good alternative methodology for ensemble quantum computing. This 
technque has many advantages: such as large working space (Floquet) that is adjustable. 
This makes a vast difference between this technique and previously proposed NMR-related 
methods because the dimension of the computing space can be augmented significantly for 
the spin systems which otherwise can only offer very low number of qubits. The sensitivity 
is generally higher than its liquid state counterpart of the same size because of the difference 
in the number densities of nuclear spins. The pulse sequences are conveniently implemented 
with conventional solid state NMR techniques although non- conventional techniques may 
be used to enhance sensitivity and achieve higher-qubit operations. We point out that the 
number of qubits currently manageable in our initial experiment is rather limited mainly 
because of the brute force pulse sequences used. We believe, however, more efficient pulse 
sequences for implementing higher-qubit QC operations in Floquet space can be found. 

The chemical shift interaction of single spin- 1/2 system is exemplified in this work, but 
the principles and procedures used here can be extended to other interactions such as dipolar 
and quadrupolar interactions, which is the object of further work. In fact, when more spins 
are involved in coupled systems or large quadrupolar interaction is present, high qubits are 
more conveniently realized in those systems. Nevertheless, we recognize that the difficulty 
involved in manipulating experimentally these systems may be significantly greater. 

Finally, for a practical quantum computer, error correction is necessary. The recently 
published error correction schemes must be modified to accommodate solid state NMR based 
quantum computing in Floquet space. 
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Appendix A: Derivation of Eq.(40) and Eq.(41) 

Prom Eqs. (24,39), the signal S^'^{t) can be written as 

S^y^^it) = ^Ti:{Ek,k'e-'''^^'Ucsk-m{\pk > +\p - e^k >)e"''^'~\< pk'\ + \p- < e^k']) (1) 

Uak'-m^iJl>Iy<l+j\e'^''^'} 
= ^Tr{E,,,,,,-,e-^('=-'='-^)-'^*C/c5fe-.„(bA; > +\p - e^k >)(< pk'\+ <p- e^k']) 

Udk'-J^>Iy<i + j\} 

In terms of matrix elements, it becomes 

Sr(t) = ^S,,,,,,S.,ne-^(^-^'-^)-*(t/^^5.-.. + U'Js7^J < n\k > 

^,AUcsk'-m~' + UlSt-m') < n'W X n'\l >< q\Iy\s ><l+j\n> 

— 2^k,k'^sA^CSk-m + ^CSk-m\Wc^Sk'-rn ^CSk'-m 

— 2^k,k' ^s,q['~^ CSk-m'-' CSk'-m CSk-m'^ CSk'-m CSk-m'-' CSk'-m '^CSk-m'^C 

Note that Ucsn is diagonal and ly is Hermitian and its diagonal elements are all zero. The 
above equation only has two non-zero terms given as 

cpm(j.\ _ ^sp \ttPP Tj-p-ep p-ep-'^ jp-e p , rfp-ep p-epj-j-pp -1 jp p^^-, /„s 

'^y V') — 2^k,k'['JcSk-ra'~^CSk'-m ^y '^'^CSk'-m '^CSk-m ^y J 

Using the matrix expression of the Ucsn{t) Eq.(34) and the orthogonal properties of Bcssel 
functions contained in coefficients A^i the time-domain signal Eq.(40) is readily obtained. 
The calculation of S'^{t) is completely the same as given above and its result is 

The quadrature detection signal (eq.(42)) is the sum of above equation and Eq.(40). 
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Appendix B: Derivation of Eq.(52) 



The effect of a gradient magnetic field is used to select certain orders of coherences from 
the initial density matrix by spatially averaging out others. An arbitrary element of Pq is 
denoted as \pi >< qj\ with coherence order of p — q and i — j in Hilbert and mode space, 
respectively. To see how gradient field select density matrix elements, we first calculate the 
following term (see Eq.(51)) 

UGjt2:tMoMjtGl:to)\pi > \qj\ (5) 

which can be written as according to Eqs.( 10,28) 

^n,m,l,kUG2n-mUioo^UG,l-k\n >< m\l >< k\pi\ < qj\e<ri-m)u.rt (g) 
— y; TI^ p;-i[uiia^,\n'>I:^<n'\+n'uirl\n'><n'\]tpTT .InTn /tjIp*^'*"'")'^'-* 

= ^n,n',mUG,.UiAr^ >< n'\e'^'''^'^UG,^-r\pm X gj|e^(»— K* 

= ^„,mUG,n-mUl^UG,m-i\pm >< qj\e^r{mt,+(n-m)t 

where tp is the 90 degree pulse width and t — t2 — to- Assuming p^{to) — \pm >< qj\, 
Eq.(51) is then written as 

^n,m,n'm'UG„Ui^UG,m-i\pm >< gmlC/^^l^^.^.t/.^lt/j^V-m'e'"'"'^"^""^'^*^^^"""'''"^'""^^*' (7) 

Let to — 0. Using Eq.(50), the terms related to the gradient field in above equation are of 
the form 

^±ilL^Guzu;rt ^±i!^Guzu>rt ^g-j 

^±iIL^G2zZU>rt ^±i^GuZUJrt 

Note that the gradient fields are symmetrical with respect to z axis and the indices 
n,m,n',m' run from — oo to oo. Therefore, above equation would vanish unless 
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(i - j)GitG, ^ (n - n')G2tG, (9) 

If we denote the coherence orders in the mode space during the first and second gradient 
field pulses, i— j and n — n', as I and k, respectively, above condition is the same as Eq.(52) 
by adding the coherence transfer condition in Hilbert space. 
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FIGURES 



Fig. 1 (a) The Floquet energy levels of spin-1/2 under MAS. (b) The read pulse (90°) for 
a Floquet state. 

Fig. 2 Upper: the sub-manifold of six energy levels of a spin-1/2 system under MAS 
and lower: the representation of each state. 

Fig. 3 The readouts of the Floquet system shown in Fig. 2 for a single crystal sample. 
The parameters are 5<75=20 kHz, rics =0.5, the spinning speed Ur—^ kHz. The relative 
orientation of the crystal with respect to the magnetic field is described by two Euler angles 
(between the lab frame and the principal axis system of the CSA tensor) (a, /3)=(30°, 60°). 

Fig. 4 The readouts of the Floquet system shown in Fig. 2 for real polycrystalline powder 
sample. The parameters: 5cs—'^^ kHz, 7705=0.5 and spinning speed ujr—^ kHz. 

Fig. 5 The pulse sequence of a 2D PASS experiment for spin-1/2 systems proposed by 
Atzutkin, Shekar and Levitt [ASL] . 

Fig. 6 The typical pulse arrangement for state labehng with gradient field selection. 

Fig. 7 Grover search based on MAS NMR. (a) Energy level diagram showing the states 
used for the experiment, (b) Schematic representation of the four Floquet states chosen for 
the experiment. 

Fig. 8 The experimental result of four possible states with Grover search algorithm 
based on MAS NMR, represented with the sideband pattern of the aromatic carbons of 
hexamethylbenzene. The rightmost peak comes from the methyl carbons and is used as 
phasing reference. 
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